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Abstract 

We investigate a zero-dimensional toy model originally introduced by Mueller and Salam [1] which 
mimics high-energy scattering in QCD in the presence of both gluon saturation and gluon number fluc- 
tuations, and hence of Pomeron loops. Unlike other toy models of the reaction-diffusion type, the model 
studied in this paper is consistent with boost invariance and, related to that, it exhibits a mechanism for 
particle saturation close to that of the JIMWLK equation in QCD, namely the saturation of the emission 
rate due to high-density effects. Within this model, we establish the dominant high-energy behaviour 
of the S'-matrix element (S") for the scattering between a target obtained by evolving one particle and 
a projectile made with exactly n particles. Remarkably, we find that all such matrix elements approach 
the black disk limit = at high rapidity Y, with the same exponential law: (S") ~ exp(— y) for all 
values of n. This is so because the S'-matrix is dominated by rare target configurations which involve 
only few particles. We also find that the bulk distribution for a saturated system is of the Poisson type. 
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1 Introduction 

Much of the recent progress in the field of high-energy QCD has come from the gradual 
understanding [2-26] of the analogies between the gluon evolution in QCD at high energy and 
a classical stochastic process, similar to the 'reaction-diffusion process' A A A [27,28]. 
Such analogies lead to the conclusion that the properties of the QCD scattering amplitudes 
at very high energy, including in the vicinity of the unitarity limit, are strongly influenced by 
gluon-number fluctuations in the dilute regime, and hence cannot be reliably computed from 
mean field approximations like the Balitsky-Kovchegov (BK) equation [5,6]. Although at a first 
sight surprising — since the high-energy regime is characterized by high gluon occupancy, and 
therefore should be less affected by fluctuations — , such a strong sensitivity to fluctuations was 
in fact noticed in early studies of unitarization in the context of the dipole picture [2-4] and, 
more recently, it has been rediscovered within the context of the non-linear QCD evolution in 
the vicinity of the saturation line [12, 14]. 

It has been realized [17] that the relevant fluctuations are missed by the JIMWLK equa- 
tion [7-9], which describes the non-linear evolution of the gluon distribution in the (high- 
density) target, as well as by its 'dual' counterpart, the Balitsky hierarchy [5] for the scattering 
amplitudes, in which the evolution is rather encoded in the wavefunction of the (dilute) pro- 
jectile. The JIMWLK equation properly encompasses the mechanism responsible for gluon 
saturation at high density, but it misses the correlations associated with gluon splitting in the di- 
lute regime. Correspondingly, the Balitsky equations faithfully describe the splitting processes 
within the projectile, but completely ignore the non-linear effects responsible for saturation. 

At the same time, new equations have been proposed [17-20], which heuristically combine 
the dipole picture in the dilute regime with the JIMWLK evolution at high density, thus leading 
to a generalization of the Balitsky hierarchy which encompasses both saturation and fluctu- 
ation effects in the limit where the number of colors Nc is large, A^c ^ 1- These equations 
have been interpreted [23] as an effective theory for BFKL 'pomerons', in which the pomerons 
are allowed to dissociate and recombine with each other, like the 'molecules' in the reaction- 
diffusion problem. (See also Ref. [29] for a related approach in the context of nucleus-nucleus 
scattering.) But the complexity of these 'Pomeron loop' equations has so far hindered any sys- 
tematic approach towards their solutions, including via numerical methods. The only properties 
of these solutions to be presently known come essentially from the correspondence with statis- 
tical physics [13-15, 17, 30, 31], which is however limited to asymptotically high energies and 
very small values of the coupling constant. 

In view of the complexity of the QCD problem, several authors [17, 32-35] have recently 
started investigating simpler models in zero transverse dimensions that allow for more direct 
studies (often analytic) of the approach towards saturation and unitarization with increasing 
energy. These models are usually formulated in terms of stochastic distribution of classical, 
point-like, particles which evolve in 'time' (rapidity) according to some suitable master equa- 
tion. The models discussed in Refs. [17, 32, 34, 35] are all borrowed from statistical physics and 
describe a. A ^ AA reaction process. The model briefly discussed in Ref. [33] has been origi- 
nally proposed by Mueller and Salam [1] as a toy model for 'dipole' evolution in the presence 
of saturation, and seems to be the closest to the actual QCD dynamics, as we shall explain at 
length in what follows. This is the model that we shall focus on in this paper. 

At this point, one may legitimately wonder about the usefulness of any zero-dimensional 
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model in the context of the high-energy evolution in QCD. As well known, this evolution is 
genuinely non-local in the transverse plane, and this non-locality is responsible for the main 
characteristic of a hadron wavefunction at sufficiently high energy, namely, the coexistence of 
two different phases at different values of the gluon transverse momentum k±: (i) a dilute 
phase at large k±, where the standard, linear (or 'leading-twist'), perturbative evolution applies 
and (ii) a high-density phase, also known as the color glass condensate [9,36,37], at rela- 
tively low k_i_, where the dynamics is fully non-linear and the gluon occupation factor (almost) 
saturates. With increasing energy, the CGC phase extends towards larger values of k±, and the 
main physical questions refer to the rate of this expansion and to the properties of the gluon 
distribution and of the scattering amplitudes in the vicinity of the saturation line — i.e., in the 
transition region from dense to dilute [38^0]. In particular, it is precisely in this region that the 
effects of the gluon-number fluctuations are most striking [14, 15, 17] : they considerably slow 
down the energy increase of the saturation momentum and render the saturation line diffuse, 
which in turn has important consequences for the measured cross-sections [41,42]. 

Clearly, all these interesting features are lost when one restricts oneself to a zero-dimensional 
problem. But even in that case, some non-trivial questions remain, whose detailed understand- 
ing could shed light on the physics of saturation in the presence of particle-number fluctuations. 
Chiefly among these, there is the problem of the approach of the (average) ^-matrix element 
(S) towards the black disk limit (S) — with increasing energy. The arguments in Refs. [1,12] 
as well as the analysis to be presented in this paper demonstrate that this approach is strongly 
influenced by the fluctuations (at least in specific frames, since the physically relevant config- 
urations are frame-dependent): the average S'-matrix near the unitarity limit is dominated by 
atypical configurations, which have a relatively low gluon occupancy — well below the average 
occupancy at that energy — and yield a large contribution ~ 1 to the average quantity (S). 

The dominance of rare fluctuations on the average -S'-matrix has an interesting implication ^ 
for the matrix element (S"") of a projectile made with n particles. As well known, the quantity 
I (5*") 1^ measures the probability that the system of n particles emerge unscattered. Simple argu- 
ments of the mean-field type suggest that (S^) ~ (S)^ at high energy. (For instance, this is the 
prediction of the Balitsky-JIMWLK equations.) However, if the average -S'-matrix is dominated 
by rare configurations for which 5' ~ 1 per incoming particle, then for these relevant config- 
urations we have 5" ~ 5 ~ 1 for any n, which after averaging yields (S^) ~ {S), in sharp 
contrast with the mean field expectation. These arguments will be confirmed by our subsequent 
analysis of the Mueller-Salam model; for instance, in the case where the target is obtained by 
evolving a single particle, we shall find that (5'") ~ exp(— F) for any n. 

Let us now explain why, in our opinion, the model introduced by Mueller and Salam in 
Ref. [1] is indeed closer to the actual QCD dynamics than the (zero-dimensional) reaction 
model A ^ A A (although the latter appears to more widely studied in the recent literature). 
The main reason is that the Mueller-Salam model correctly captures the QCD mechanism for 
the saturation of the particle number at high energy, which is not particle recombination (as in 
the reaction-diffusion process), but rather the saturation of the rate for particle emission due to 
high-density effects. The new particle which is emitted when increasing the rapidity in one step 
can be radiated off any of the particles created in the previous steps and, moreover, it can un- 
dergo multiple scattering off the latter. When the density is relatively low (the low-energy case). 



^ To our knowledge, this implication has not been previously noticed in the literature. 
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the emission rate is simply proportional to the total number of sources, leading to an exponen- 
tial increase of the number of particles with the rapidity Y. But when the density becomes large 
enough for multiple scattering to be important, the emission rate saturates at a constant value, 
independent of the number of particles in the system. Then, the number of particles keeps grow- 
ing, but only linearly in Y. This 'almost saturation' scenario is indeed similar to the physical 
mechanism at work in QCD [43,44], where gluon saturation — as described by the JIMWLK 
equation — proceeds via the saturation of the rate for gluon emission by strong color field ef- 
fects (see also the discussion in Refs. [10,45,46]). By contrast, in the reaction-diffusion process, 
the particle (occupation) number saturates at a fixed value (independent of Y), a situation which 
looks unphysical from the perspective of QCD: when increasing the rapidity even within the 
saturation domain, the gluon phase-space opens towards lower longitudinal momenta, which 
then allows for further radiation. 

Closely related to the above, there is a second reason why the Mueller-Salam model is better 
suited to mimic the QCD evolution than the reaction model : unlike the latter, the former is 
consistent with the boost invariance of the S'-matrix in the presence of multiple scattering. 
In fact, in Ref. [1], this model has been constructed precisely by requiring that the (average) 
S'-matrix be independent upon the choice of a Lorentz frame. Remarkably, within the zero- 
dimensional context, the requirement of boost invariance together with the assumption of one 
particle emission per unit rapidity are in fact sufficient to uniquely fix the rate at which a system 
of n particles can emit another one, and predict, in particular, the saturation of this rate for 
sufficiently large values of n. On the other hand, a model including both particle splitting {A 
AA) and particle recombination {AA A), is in conflict with boost invariance as we shall 
explain in Appendix A. 

In turn, these conceptual differences between the Mueller-Salam model and the reaction- 
diffusion model entail important structural differences between the two, as well as differences 
between their respective predictions. For instance, we shall see that the evolution described by 
the Mueller-Salam model involves particle-number-changing vertices m ^ n for all values of 
m and n, so like the JIMWLK equation and its recent generalizations in Refs. [21,22,24,25]. 
By contrast, the reaction model involves only 1 — > 2 and 2 — > 1 vertices. This difference is 
in fact related to the issue of boost invariance in the presence of multiple (eikonal) scattering: 
As recently shown in Refs. [22,23], the requirement of boost invariance implies a symmetry 
property for the evolution Hamiltonian known as self-duality. As further explained in Ref. 
[23], on the example of the 'Pomeron loop' equations, the reaction-diffusion dynamics can 
be made consistent with self-duality, but only under the additional assumption that one can 
ignore multiple scattering for the individual particles — an assumption which looks unnatural 
in the context of high-energy QCD. By contrast, if one allows for multiple scattering in the 
eikonal approximation, then a self-dual evolution Hamiltonian must involve an infinite number 
of gluon-number-changing vertices, as explicit in the QCD constructions in Refs. [24,25]. 

Note also another feature of the reaction-diffusion models (at least in zero transverse di- 
mensions), which looks unphysical in the context of QCD: As shown in Ref. [35], such models 
predict a 'grey disk' limit at high energy, that is, the corresponding S'-matrix saturates at a 
non-zero value when Y oo. This will be briefly discussed in Appendix A. 

Other interesting results emerging from our study, may have a counterpart in QCD as well. 
The nature of the physically relevant configurations which contribute to the average S'-matrix 
is frame-dependent. For definiteness, consider the collision between two systems which in their 
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respective rest frames reduce to one particle (the 'target') and, respectively, to n particles (the 
'projectile'). While the average ^-matrix (S*") is frame-independent, the configurations which 
dominate the average at high energy are different in different frames: whereas in the rest frame 
of the projectile, the average is dominated by rare target configurations which are dilute, in the 
target rest frame, and for n > 2, {S'^) is rather controlled by the typical configurations in the 
wavefunction of the projectile, which are at saturation. We shall find that these typical config- 
urations at saturation follow a Poisson distribution. The case n = 1 (the symmetric scattering) 
is special^ : (S) is controlled by rare configurations, but marginally sensitive to the saturated 
ones, in any frame, and has a different subleading behaviour at high energy as compared to the 
case where n > 2. 

The paper is organized as follows: Sect. 2 will be devoted to structural aspects of the Mueller- 
Salam model. In Sect. 2.1, we shall construct the master equation from the requirement of 
boost invariance; then, in Sect. 2.2, we shall use this equation to derive a hierarchy of evolution 
equations for the yS-matrix elements (5*"). In Sect. 3 we shall use the solution to the master 
equation to investigate the high-energy behaviour of the model. We shall successively consider 
the S'-matrix elements {S^) (in Sect. 3.1), the bulk of the particle distribution (in Sect. 3.2), 
and the frame-dependence of the configurations which control the average S'-matrix (in Sect. 
3.3). Finally, in Sect. 4, we shall discuss the correspondence between the toy model and the 
high-energy problem in QCD, with emphasis on the similarity between the respective evolution 
equations in various limits. 

2 The toy model: Structural aspects 

As explained in the Introduction, the toy model that we shall consider is 'zero-dimensional' 
in the sense that its only variable is the rapidity interval Y which controls the high-energy 
evolution. The 'physical' problem that we have in mind is that of the scattering between two 
systems of particles, referred to as the projectile and the target. Each system is characterized by 
a probability distribution giving the number of particles it contains at a given rapidity. During 
the collision of the two systems, particles of the projectile and the target undergo independent 
collisions, characterized by a scattering amplitude r. In order to keep close contact with QCD, 
we shall use throughout a QCD-inspired terminology and refer to the two colliding systems 
as two "onia" and to the particles which compose these onia as "dipoles". This is suggestive 
since, in the dilute regime at least, the evolution that we shall describe corresponds to the zero- 
dimensional version of Mueller's dipole picture [1^]. 

The number of dipoles in each system depends upon its rapidity, and therefore on the frame. 
We put the right mover (the 'target') at rapidity Y — Yq and the left mover (the 'projectile') at 
rapidity — Yq (^o > 0). We denote by P^^{Yq) the probability to find exactly m dipoles inside 
the projectile at rapidity Yq. Similarly P^iY — Yq) will be the corresponding distribution for 
the target. Allowing the dipoles in each onia to scatter independently with each other, the S'- 
matrix for a configuration with m dipoles in the projectile and n dipoles in the target takes the 
factorized form S — cr"*", where a — 1 — r is the ^-matrix for the scattering between two 



Our result for (5) will be the same as in the original analysis in Ref. [1], where however the precise 
nature of the relevant configurations has not been investigated. 
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elementary dipoles and r is the corresponding T-matrix. The most interesting case in view of 
the comparison with QCD is the weak coupling regime r ^ 1 (see Sect. 4). 

The physical ^-matrix is obtained by averaging over all the possible configurations in the 
two onia, with the respective probability distributions: 

oo 

{S)y = E P^iYo)P^{Y - Yo) <7-" (2.1) 

m,n=l 

This expression for the S-matrix reflects the fundamental factorization assumption which lies at 
the basis of our analysis. This formula was used in Ref. [33], while in their original formulation, 
Mueller and Salam [1] write, in their analog of Eq. (2.1), e"™*^ in place of cr™". 



2.1 Particle saturation from boost invariance 



We shall assume that the two onia follow the same evolution law in rapidity, that is, they 
obey the same microscopic dynamics. The initial conditions at F = are left arbitrary for the 
time being, but they will be specified later, when this will be needed. Then, as we show now, 
the rapidity evolution is uniquely fixed by the following two constraints: 

(i) Lorentz (boost) invariance: The total yS-matrix should be independent of the choice of 
frame, i.e. , of Yq, which implies 



= 0. (2.2) 



(ii) The onium evolves through dipole emission, in such a way that only one new dipole can 
be produced under a step dY in rapidity. Thus, in a step dF, a system of n dipoles can turn 
into a system of n + 1 dipoles, with a probability /„dy, and stay in its initial configuration 
with a probability 1 — /ndF. It follows that the probability -Pn(^) evolves according to the 
following master equation: 

= P„_i(F) - fnPniY), (2.3) 

where /„ is a function of n to be determined shortly. In this equation and throughout, Pn 
refers generically to either P^, or P^. Note that the evolution (2.3) preserves the total 
probability: d(En Pn{Y))/dY =0. 

As noticed in Ref. [1], the two conditions above determine the transition probabilities /„ up 
to an overall normalization factor (which is then fixed by /i). Indeed, setting the derivative of 
Eq. (2. 1) with respect to Yq equal to zero and using Eq. (2.3) we arrive at 

P^(Yo)P^{Y - Yo) - <7") - Ul - a-)] = ^ = c (1 - a"), (2.4) 

m,n 

where c is a constant that can be reexpressed as c = — c), where fi is the rate of splitting 
of one dipole into two. Without loss of generality, we shall set fi = 1. (Within the QCD dipole 
picture, this rate would be equal to (Xsj^ix times the dipole kernel; see, e.g., Refs. [2, 11] for 
details.) We then obtain, from Eq. (2.4), 
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1 — cr" 

fn = 1 = l + a + --- + (7"-^ < n. (2.5) 

1 — a 

The quantity /„ is the rate at which a system of n dipoles emits an extra one in a step dY in 
rapidity. The constraint of boost-invariance of the scattering matrix forces the /„ (and hence 
the probability distribution P„) to depend on a, and furthermore induces non trivial correlations 
that affect their n-dependence. To better visualize the implications of Eq. (2.5), it is useful to 
consider the two limits of strong (r — > 1) and weak (r 0) couplings. In the weak coupling 
regime, cr ^ 1 and ~ n: in this regime, the n dipoles split independently of each other. In the 
strong coupling regime however, a ~ 0, and /„ ~ 1: in this case correlations in the n-dipole 
system are such that only a single dipole can be emitted. 

When T is fixed and small, which is the case of physical interest, a similar change of regime 
occurs when increasing the value n. Namely, Eq. (2.5) implies that, when r <^ 1, the rate /„ 
has the following limiting behaviors: 



fn 

n 



^_T(ril)^___ forn«l/r 

1 2 (2.6) 
• • • < 1 for n > l/r. 



rn 



According to Eq. (2.6), the emission rate fn saturates at a large, but constant, value 1/r 3> 1 
when n ^ rigat ~ l/r ^ 1. As we shall discover soon, Usat is indeed the order of magnitude of 
the average dipole number at the onset of saturation. 

Alternatively, observing that cr" is the yS-matrix for the scattering of a dipole projectile on 
an n-dipole target, we may interpret the presence of higher powers of a in the emission rate 
(2.5) as reflecting the multiple scattering between the newly emitted dipole and its sources. This 
interpretation becomes perhaps more transparent when Eqs. (2.3) and (2.5) are used to compute 
the rate for the evolution of the average ^'-matrix with Y. One finds 

^ = -- E Pl:.{yo)P^{Y - Y,) a-"(l - a-)(l - a"), (2.7) 

m,n 

where the quantities t„ = 1 — c" and t„i = 1 — o-™ can be recognized as the scattering amplitudes 
for the scattering between one dipole — the one created in the last step in the evolution — and 
the n preexisting dipoles in the target and, respectively, the m dipoles in the projectile. 

For latter reference, let us notice that Eq. (2.1) is equivalent to the following factorized form 
of the scattering amplitude (T) = 1 — (S) 

{T)y = E S:^(^^'Vo {t')Y-Yo, (2.8) 
k=l 

which is the form generally used in the QCD context (see, e.g., Refs. [6, 16,26]). In this equation, 
(n*^*^)) = (n(n — 1) • • • (n — /c + 1)) is the /c-body normal-ordered dipole number (here, in the 
wavefunction of the projectile), t = 1 — s the scattering amplitude for a single dipole, and 
{t^)Y-Yo is the average amplitude for the simultaneous scattering of k dipoles off the target. 
Note that, in the discussion above, we have introduced the notation s for the 5'-matrix for a 
projectile made with a single dipole. This notation will be systematically used in what follows. 
Correspondingly, the S'-matrix for a projectile made with exactly m dipoles is s'^. 
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Finally, let us mention that the P„'s can be obtained from the following generating functional 

oo 

Z{u,Y)^Y.Pn{Y)u''. (2.9) 



n=l 



from which most quantities of interest can also be derived. In particular, 

d''Z{u,Y) 



(n^'^)) = (n(n - 1) • • • (n - A; + 1)) 



(2.10) 

u=l 



Differentiating Eq. (2.9) with respect to Y and using the master equation one arrives at 

^ = \^[Z{au)-Z{u)]. (2.11) 

This equation turns out to be difficult to solve analytically, so in what follows we shall work 
mostly with the master equation. 

2.2 Evolution equations: Saturation, unitarity & fluctuations 

Using the master equation (2.3), we shall now deduce evolution equations for physical quan- 
tities such as the average S-matrix or the average dipole number, and study some of their gen- 
eral properties. More specific predictions about the solutions to these equations at high energy 
will be discussed in Sect. 3. 

Let us begin with the scattering matrix and rewrite Eq. (2.1) as 

(^>y = E^m(^o) (Oy-yo (2-12) 

m 

where 

{s'^)y = J2 PniY) f^™" = ^'J",Y) (2. 13) 

n 

is the average S'-matrix for a projectile made with exactly m dipoles and a generic target. Matrix 
elements like (s"*)y carry information about the dipole correlations in the target. 

From Eqs. (2.13), (2.3) and (2.5), it is now straightforward to obtain the following evolution 
equation for (s"*) : 

-^ = fm[{s"'^')-{sn]- (2.14) 

This is not a closed equation — (s"^) being related to — , but rather a particular equation 

from an infinite hierarchy. This equation has been obtained here by following the evolution of 

the target (namely, by using Eq. (2.3) for P^), but it can be easily reinterpreted as describing 
evolution in the projectile: when increasing the projectile rapidity by dY, the incoming system 
of m dipoles can turn, with a rate f„i into a system of (m + 1) dipoles, which then scatters off 
the target with an ^-matrix (5™+^). This is the origin of the first term within the brackets in the 
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r.h.s. of Eq. (2.14). As for the second term, involving (— (s"^)), it corresponds to the case where 
the system remains intact during the evolution (which occurs with probability 1 — fmdY). 
For m = 1 (a single dipole in the projectile), Eq. (2.14) reduces to 



dY 



is') -is), 



(2.15) 



which is formally identical to the first equation in the Balitsky-JIMWLK hierarchy. However, 
differences with respect to the latter appear already for m = 2 . Indeed we have 



dY 



= (2 - r) is') - {s') 



(2.16) 



The corresponding Balitsky equation^ would not contain the term proportional to r in the 
r.h.s. More generally the Balitsky equations are obtained by replacing by m in the r.h.s. of 
Eq. (2.14), which amounts to ignoring saturation effects in the projectile. Such effects may be 
expected to be negligible when mr <^ 1, but we shall discover in the next section that this is 
not so. In fact they are essential to get the correct description of the evolution at large Y. 

Alternatively, the difference between the hierarchy in Eq. (2.14) and the Balitsky hierarchy 
can be attributed to particle-number fluctuations in the target. For n 1 the discrete nature of 
n becomes inessential and the associated fluctuations are unimportant. Therefore, let us replace 
the summation over n in Eq. (2.1) for the average <S'-matrix by the corresponding integration 
(but keep a discrete sum over m on the side of the projectile) and require boost invariance. 
Then we find that we need to impose a separate evolution law for the target and the projectile 
wavefunctions, that is, the splitting rates and must be different functions. It is a matter of 
simple algebra to obtain ^ 



dY 



d_ 
dn 



a 



Incrl 



d_ 
dn 



fnPniy) 



(2.17) 



while at the same time the left mover is evolving according to Eq. (2.3), but with = m. After 
also trading the sum over n in the definition (2.13) of (s™) by the corresponding integral, we 
find that Eq. (2.17) leads to the Balitsky hierarchy, as anticipated: 



d{s'' 



dY 



m 



(2.18) 



The fact that neglecting the correlations associated with particle-number fluctuations in the 
target is equivalent to ignoring the saturation effects in the projectile is in agreement with the 
general arguments in Ref. [17]. Further insight on this issue can be gained by rewriting the 

evolution equations in terms of the scattering amplitudes (t^) introduced in Eq. (2.8). The cor- 
responding equations are easily deduced from those for (s"*) by using t = \ — s. The first three 
equations in this hierarchy read: 



^ In the remaining part of this section, by the "Balitsky equations" we shall mean the simplified version 
of these equations at large- A^c and for zero transverse dimensions, that is, Eq. (2.14) with m. 
^ In Sect. 4, we shall argue that Eq. (2. 17) is the toy-model analog of the JIMWLK equation. 
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d{t) 



it) 



it'), 



(2.19) 



dY 



2{t')+T{t{l-tf), 



(2.20) 



dY 



d{t') 



?,{t^) + ?>T(f{i-e)) + T^{t{i-tf). 



(2.21) 



dY 



The new terms, as compared to the corresponding Balitsky equations, are those proportional 
to r or in the last two equations. As mentioned before, these terms reflect dipole-number 
fluctuations in the target. 

First we consider Eq. (2.20): Among the three fluctuation terms there, namely r(t(l — t)^) = 
T{t^ — 2t^ + t), we shall focus on the last one, r(t), since this is the most important one^ . 
(When T <^ 1, the other terms are negligible in all regimes.) Although formally suppressed by a 
power of r with respect to the (BFKL-like) term {t^), the fluctuation term r(t) is in fact equally 
important in the low density regime where (t) ~ r. To clarify its physical meaning, note that, in 
the dilute regime, the average dipole scattering amplitude is simply proportional to the average 
dipole number in the target: (t) ~ r (n). Thus, r(t) ~ T^(n), and the physical interpretation 
becomes transparent: under a rapidity step dY, any one among the (n) target dipoles can split 
into two and then the child dipoles can scatter with the two projectile ones, with strength r^. 
Thus, the simultaneous scattering of two projectile dipoles gives us access to the correlations 
induced via dipole splitting in the target. 

Let us now move to the next equations in the hierarchy. Eq. (2.21) for (t^) contains two 
relevant fluctuation terms, of order r(t^) and T^(t), respectively. The first one has the same 
physical origin as the term r(t) in the equation for {t^): that is, two among the three projectile 
dipole 'feel' a fluctuation by scattering off the child dipoles produced by a splitting in the 
target. The second term, r^(t) ~ T^{n), describes the process in which the fluctuation is felt 
by all the three projectile dipoles: e.g., two of them scatter off one child dipole, and the third 
one scatters off the other child dipole. These terms are both relevant since they are of the same 
order in the low density region where (t) ~ r. It is not hard to understand the generalization to 
higher equations in this hierarchy: the equation for (t*^) will involve relevant fluctuation terms 
of the following types: r(t^~^), r^(t^~^), . . . , T^~^{t). All such terms appear to be important 
for building up many-body correlations in the dilute regime. 

Incidentally, the previous discussion also shows that the r.h.s. of the evolution equation for 
{t^) involves terms proportional to {f^) where the power m can take all the values from m — 1 
up to m = A; + 1. Since k is generic, this means that the hierarchy as a whole involves vertices 
relating {t'^) to (t™) for arbitrary values of k and m. 

Eq. (2.19) allows us to estimate the rapidity Yc for the onset of unitarity corrections in the 
dipole-target scattering. Let us assume that, at y = 0, the target consists in a single dipole: 
P^(0) = 5ni, then we have {t^)^ — t^, hence {t^'^^)^ <S {t^)^ and these inequalities will be 
preserved in the early stages of the evolution where one can therefore neglect (t^) as compared 
to (t) in the r.h.s. of Eq. (2.19). This then leads to (t) ~ re^, the analog of the BFKL increase 
[47]. Clearly this ceases to be correct when (t) ~ 1, that is, for Y ^ Yc = ln(l/T). For larger 
rapidities Y > Yc, multiple scattering becomes important and ensures the unitarization of the 

^ This is the analog of the 'fluctuation terms' which appear in the Pomeron loop equations constructed 
in Refs. [17, 18]. See the discussion in Sect. 4. 
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dipole amplitude ((t) 1 as F — > cxd), as we shall discover in the next section. 

Equivalently, marks the onset of the saturation effects in the target in the frame where the 
projectile is dilute. To see this, consider the evolution equation for the average dipole number 
(n) in the target, that is 

= - (1 - a ) , (2.22) 

where (1 — a") = (t) is recognized^ as the scattering amplitude for a projectile made with 
one dipole, cf. Eq. (2.13). At low density, (1 — a") ~ t (n), and the dipole number exhibits 
the exponential growth (n) = characteristic of BFKL evolution. But when (n) becomes as 
large as Ugat ~ I/t, which happens for Y ~ Yc, the growth is tamed by non-linear effects. 
Eventually, when Y » Y^ (cr") becomes negligible so that d(n)/dy ~ 1/r and (n) grows 
only linearly (as expected for the gluon occupation factor in QCD [43,44]). To summarize, 

(n) = when F < Fc ; (n) ~ when F > Fc ■ (2.23) 



3 The high-energy behaviour of the toy model 



In this section, we shall establish the dominant high-energy (i.e., large-F) behaviour pre- 
dicted by the toy model for the dipolar ^'-matrix elements (s"^) with m > 1, and also for the 
dipole distribution in the target. 



3.1 The dipole S-matrix elements {s"^) 



For definiteness, let us focus on the situation where the target reduces to a single dipole in its 
own rest frame: -P^(O) = Sni- (More general initial conditions will be briefly discussed in Sect. 
3.3.) Also, let us perform our calculation in the projectile rest frame. This last choice turns out 
to be very non-trivial since, as we shall discover in Sect. 3.3, the configurations which dominate 
the average ^'-matrix are different in different frames. In conformity with these choices, the 5*- 
matrix element (s™) will be computed according to Eq. (2.13); that is, we shall first solve the 
master equation (2.3) with the initial condition -P^(O) = Sni and then perform the sum over n 
in Eq. (2.13). We shall not be able to perform these operations exactly, but only under suitable 
approximations, which preserve the dominant behaviour of (s™) at high energy. 

Before we proceed with our calculations, we note that the hierarchy in Eq. (2. 14) admits the 
following, two-parameter, family of solutions: 

{s).. = - (F-F„)e-^ and (Oas = —e'"^ for m>2, (3.1) 

1 — 0" 1 — (J™- ^ 

with arbitrary values for the parameters A and F^. Although these cannot be the complete so- 
lutions, as they do not obey the physical initial conditions at F = for any choice of the free 
parameters, with the choices A = a"^ and F^ ~ Fc = In (1/r) they describe the asymptotic form 
of the physical solutions at large F » Fc (hence the subscript 'as'), as we shall see. 

^ Alternatively, the r.h.s. of Eq. (2.23) is recognized as the average emission rate (/„). 
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What is remarkable about the behaviour of these solutions is the fact that, for asymptotically 
large Y, all the 's approach the black disk limit (S = 0) according to the same exponential 
law exp(— F), for all m. This might look counterintuitive since, at large Y, the target onium is 
typically characterized by a large number of dipoles, off which a projectile dipole will scatter 
with a very small S'-matrix, s <S 1. We would then expect (s"*) <^ (s) when m > 1. In fact, if 
one assumes that {s"^~^^) <^ {s"^) for sufficiently large Y, then one gets from Eq. (2.14) 

(s")typicai ~ e-^-^ for large Y , (3.2) 

which is indeed consistent^ with the initial assumption <g (s"*) (at least, as long as 

m <^ 1/t), but is nevertheless in contradiction with the correct asymptotic behaviour exhibited 
in Eq. (3.1). The problem with the estimate (3.2) is that it implicitly assumes that the average 
^-matrix is controlled by the typical target configurations, which have a large dipole number. 
However, even at large Y there is still a non-vanishing probability that dilute configurations 
remain present in the target, and as we shall show explicitly, the sum over n in Eq. (2.13) is 
dominated by rare fluctuations for which n is relatively low, n ~ 0{1), and therefore s ~ 
a" ~ 1, which in turn implies s™ ~ s. Such fluctuations have a low probability ~ exp(— F), 
but this is compensated by the fact that their contribution to the average ^-matrix is relatively 
large, of order one. On the other hand, the typical configurations have a probability of order 
one. Since n is very large and of the order of the average dipole number (n) ~ (Y — Yc)/t (see 
Eq. (2.23)), the contribution of these configurations to (s™) is exponentially suppressed, namely 
cr™" ~ exp{— m(F — K)}. These estimates suggest that indeed, at least for m > 2, the average 
^-matrix (s™) is dominated by the rare dilute configurations. The case m = 1 is a priori more 
subtle, since then both rare and typical configurations can give significant contributions. 

To make this discussion more quantitative and confirm the asymptotic behaviour displayed 
in Eq. (3.1), we shall now proceed to an explicit calculation of the probabilities Pn{Y). We do 
so by using the Laplace transform of -Pn(F) 



oo 



Pn{ou) = I dYe-^'^PniY), (3.3) 



in terms of which the master equation reads 

Using the initial condition Pn{0) — Sni, we get 

W = I fi f 1 + f 1 \ (3-5) 

Jn k=l \ JkJ 



^ But, clearly, this consistency disappears for large m ^ 1/t, which explains why the Ansatz (3.2) 

cannot be a solution to the complete hierarchy. In fact, because of the coupling between small and large 
values of m throughout the hierarchy, the estimate (3.2) is wrong even for relatively small m > 2, as 
clear from the comparison with Eq. (3.1). 
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from which one can obtain Pn{Y) by the inverse Laplace transform 

duj 

c 



Pn{Y)^ /^e-^P„(a;), (3.6) 



where the integration is to be done in the counterclockwise direction along any contour C enclos- 
ing all the poles of -P„(c^). Notice that all the these poles occur in the finite interval (— l/r, — 1]. 

There are two limiting behaviors of Pn{Y) that can be identified respectively with the cases 
of weak coupling (/„ = n) and strong coupling (/„ = 1). In the first case one obtains the 
familiar distribution of the dipole model [3] 



P^^ny) - - e-^y-\ (3.7) 

In the second case, the distribution is of the Poisson type: 

pPoisson(y)^_ZL_e-^. (3.8) 

In the general case, we do not have a closed expression but one can easily construct Pn{Y) in 
the form of a (finite) sum: 

n 

P„(r) = ^ c^-^'=^ (3.9) 

k=l 

with coefficients cl determined by 

n r 

Y.cl = 6ni and cl = -/^cl-\ (3.10) 
k=i 

These can be iteratively constructed. Here we shall present only few of them, some of which 
will be important for our subsequent analysis. We have 

^ n ^ fn-l n ^ (-l)"/n-l n ^ (-l)""^^ H in 



Just for illustration, the first four probabilities are given by 

Pi(y)=e-^, (3.12) 

P2(F) = le-^-ie-^^^, (3.13) 

Pm = A e-^ - e-^^^ + e-^^^, (3.14) 

^Y) = ^ - e-/^^ + e-/3^ - \ e"/^^. (3.15) 

It is easy to verify that these distributions go over to either the dipole distribution (3.7) or the 
Poisson distribution (3.8) as cr ^ 1 or cr ^ 0, respectively. 



2rt-3' (jn(n-l)/2 ' " ^n(n-l)/2 
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Fig. 1. For high Y and values of n such that n < ricr the probabiUties Pn{y) are dominated by the first 
term of the sum in Eq. (3.9). Notice that for the values ofY and r used in this plot, ricr ~ 67. 

From Eq. (3.9)-(3. 1 1), it is clear that for large Y and not too large values of n, the dominant 
contribution to P„{Y) is given by the first term in Eq. (3.9), proportional to e~^. Indeed, the 
terms with > 2 are exponentially suppressed at large Y with respect to the first term. But 
the situation changes when n becomes large, since the coefficients increase rapidly with n 
and this rise can compensate for the exponential suppression with Y. One can roughly estimate 
the value of n at which this change of regime occurs by requiring that the first two terms in the 
expansion of PniY) become of the same order. This criterion yields 

aY - ln(l/r) Y-Y, 

Ucr ~ n \ ~ ' (3.16) 

|m(T| r 

where Yc = ln(l/r), as before, and the second estimate holds when r ^ 1. A more pre- 
cise estimate for ricr will be obtained in the next section and reads ricr ~ (cr/r)(F — Yc). 
Clearly, this number is of the order of, but smaller than, the average number of dipoles (n) 
at large Y (cf. Eq. (2.23)). Hence, the configurations with n <^ n^r are relatively dilute and 
thus have a small probability, approximately given by the first term in the sum in Eq. (3.9) : 
PfiiX) ~ (l/o")"^^e^^ <^ 1. In Fig. 1, this analytic estimate is compared with the exact result, 
as obtained by the numerical solution to the master equation. On the other hand, for ricr — 
the case of the typical configurations — PniX) is of order one and is dominated by the terms 
with large values of k, of order n (see Sect. 3.2). As discussed above, we do not expect these 
bulk configurations to affect the calculation of the average S-matrix elements (s™), which are 
expected to be controlled by the dilute configurations with n <^ Ucr. 

Let us verify this by explicitly computing (s™), starting with m = 2. To that aim, we separate 
out the contribution of the configurations with n < Ucr in Eq. (2.13) for (s^). This yields 

{s') = V Pn{Yy^ ^ V ^ = ^AL^ le-y ^ ^e-^, (3.17) 
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Fig. 2. The distribution a'^^Pn as a function of n for large Y; only ttie rare configurations with n up to 
~ l/r contribute to (s^). 

where the neglected terms are of (9(e~^*^^~^'^)), and thus are exponentially suppressed when 
F ^ Fc relatively to the dominant contribution, of C(e~*^^~^'=^). In practice only the configu- 
rations with up to ~ l/r dipoles contribute to the final result (this is manifest in Fig. 2), and 
for them cr" ^ cr^/^ ~ e~^. For any such configuration, the S-matrix for a projectile dipole is 
of order one, s = a" ~ as anticipated. The contribution of the bulk configurations with 

n ^ Tier to (s^) will be considered in Sect. 3.2 and found to be of (9(e~^*^^~^'=''), so like the 
terms neglected in evaluating Eq. (3.17). Hence, the final result in Eq. (3.17) gives indeed the 
dominant behaviour when F ^ Y^. 

One can extend this calculation to an arbitrary m > 2. In fact, the larger m is, the faster the 
convergence of the sum over n in Eq. (2.13) is . Then it is straightforward to show that 



cr 



-Y 



a 



m—l 



for 



m > 2. 



(3.18) 



We emphasize here that not only the F-dependence, but also the prefactor in the above equation 
are well under control. This result is in agreement with the respective one in Eq. (3.1) and it 
fixes the parameter A there to be A = cr^. 

More generally, the above procedure allows one to calculate the generating functional Z(u, F), 
Eq. (2.9), for large F and any value of u which is strictly smaller than a: the corresponding re- 
sult is obtained by simply replacing cr'" u in Eq. (3.18), and reads 



Z{u,Y) 



u 



-Y 



1 



[u/a) 



for F > Fc and < m < cr . 



(3.19) 



The analytic estimate (3.18) for (s™) and the corresponding one for Z{m) 
compared to the respective numerical results in Figs. 3 and 4. 



Ziu 



^'^ The extreme limit of this behaviour occurs when m becomes larger than l/r. Then the only configu- 
ration of the target wavefunction which is relevant is the one with n = 1. 
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Fig. 3. The S-matrix for the scattering of m projectile dipoles normalized to that for the scattering of 
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Fig. 4. The generating functional as a function of m = Inu/lno" for m > 1. When m is an integer, 
Z{m) gives the S-matrix for the scattering of m projectile dipoles off the target. 

Let us now turn to the case m = 1, which is special. Then the analog of Eq. (3. 17) reads 

ncr "cr _2 

{s) = Y: Pn{Y)a- ^ ae-^ E 1 ^ 7^ - ^c) e"^ , (3.20) 

n=l n=l -'- ^ 

where we have used the improved estimate ricr ~ (cr/r)(F — Yc), to be found in Sect. 3.2. The 
prefactor in front of the exponential in the final result is essentially the number of configurations 
which contribute to the average S-matrix, with each such configuration bringing a contribution 
of 0{e^^). Note, however, that the above sum is dominated by its upper limit, i.e., by config- 
urations with n ~ ricr, for which our approximations are not fully under control. This result 
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Fig. 5. The distribution a^Pn as a function of n for large Y; all the configurations with n up to 
Ucv ~ o'{Y — Yc)/t contribute to (s). 
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Fig. 6. The expectation value {s) of the S-matrix for onium-onium scattering. Starting from {s)o = a, 
this expectation values falls exponentially to zero when Y » Yc, in agreement with the analytic predic- 
tion (3.20) (see also Fig. (7)). 

turns out to be nevertheless correct (up to corrections of 0{1) to the rapidity shift Y^, which 
go beyond the present accuracy), for the following reason: the distribution cr"P„(F) is almost 
flat as a function of n so long as n < ricr — as manifest on Eq. (3.20) — but it drops out very 
fast when n > (this can be seen in the numerical results in Fig. 5 and will be analytically 
verified in the next section). 

A different way to check Eq. (3.20) is to rely on the hierarchy of evolution equations for 
(s™). As previously noticed, (s™) given by Eq. (3.18) is an accurate asymptotic solution for 
all the equations in the hierarchy starting with the second one. Then one can use this result, 
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Fig. 7. The average value (s) for higher values ofY. The analytic curve is obtained by using Eq. (3.20) 
withY^ = ln{l/T). 

together with the first equation in the hierarchy, Eq. (2.15), in order to determine the asymptotic 
form of (s). Thus one easily recovers Eq. (3.20) where however the rapidity shift Yc is left 
undetermined. A further confirmation of Eq. (3.20) can be found in the original calculation of 
(s) by Mueller and Salam [1], which is based on the solution to Eq. (2.17) and which yields 
the same result as our Eq. (3.20) with Y^ = ln(l/r) + The analytic estimate (3.20) is 

compared to the respective, exact, numerical result in Figs. (6) and (7), which confirm that both 
the F-dependence and the normalization shown in Eq. (3.20) are indeed under control. 

To conclude this section, it is instructive to compare the above results to those predicted by 
the dipole picture. Since the latter is not boost invariant, we shall obtain different results for 
(s™) depending upon the system that we decide to evolve: the target or the projectile: 

(i) Projectile evolution in the dipole picture. The dipole ^-matrix elements (s™) obey the 
Balitsky equations, that is, Eq. (2.14) with the prefactor fm replaced by m. By doing this re- 
placement, one looses Lorentz invariance. This is formally recovered by allowing the projectile 
and the target to obey different evolution equations. Namely, if the projectile obeys dipole evo- 
lution, then in order for the r.h.s. of Eq. (2.12) to be independent of Yq, the target probabilities 
P^{Y - Yo) which enter via Eq. (2.13) must evolve according to the 'JIMWLK' 

version of the master equation , i.e., Eq. (2.17). 

Let us then estimate the high-energy behaviour of (s™) as predicted by the Balitsky equa- 
tions: by inserting the Ansatz (s™) = exp (— CmF) in this hierarchy, one immediately finds 
Cm = m, and hence 

(s"^) = e~™^ at large F (no saturation in the projectile) , (3.21) 

which is essentially the mean field estimate (3.2). Thus, by neglecting saturation effects in the 
projectile wavefunction (or, equivalently, particle-number fluctuations in the target wavefunc- 

This is, of course, in agreement with the fact that the reduced Balitsky equations (2.18) can be also 
obtained from the 'JIMWLK' evolution (2.17) of the target. 
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tion), we are led to the wrong conclusion that (s™"^^) vanishes faster than (,s™) at large Y. 

(ii) Target evolution in the dipole picture. By using the target-average expression (2.13) 
for (s™) together with the explicit solution (3.7) for the probabilities in the dipole picture, one 
immediately finds: 

{s^) ~ — e"^ for m > 1 and large Y (no saturation in the target), (3.22) 

where the neglected terms are of order (9(e~^^). For m > 2, this is essentially the same as 
the correct result at large Y , as given in Eq. (3.18). This agreement is consistent with the fact 
that (s™) with m > 2 is dominated by rare target configurations with a small dipole number 
and thus are insensitive to saturation effects. For m — 1, on the other hand, the dipole-picture 
prediction in Eq. (3.22) is different from the correct respective result in Eq. (3.20); namely, it 
misses the large, overall, factor Y — Yc, which in Eq. (3.20) has been produced by summing 
over configurations with a relatively large number of dipoles n ^ (Y — Yc)/r (cf. Eq. (3.20)), 
which are close to saturation. 

3.2 The bulk of the dipole distribution 

Although they appeared to be quasi-irrelevant in our previous calculation of the average 
S'-matrix, the typical configurations become important when this calculation is performed in 
a different frame. Therefore, in this subsection we study these typical configurations, i.e. the 
probability distribution Pn{Y) at high energy (typically, Y > Yc) and for a large number of 

dipoles n > risat- 

To compute this distribution, we return to the exact Laplace transform, Eq. (3.5), and write 

ln[/„P„(^)] = -Eln l + f =E E (3-23) 

k=l \ Jk/ r=l ' fc=l 

where we have expanded the logarithm and exchanged the order of the two summations. We 
shall now perform the summation over k in the above equation, under the assumptions that 
nr » 1 and r <^ 1. More precisely, we shall neglect terms which are of order (or smaller than) 
C(e~^") and/or 0{t). Notice that all these summations grow linearly with n for large n (since 
/^^ ~ r*" when A; ^ 1/r), so it will be convenient to subtract this large contribution and then 
perform approximations on the remainder. We thus have 

El = E fk' = ^ E TZ-n^ - ™ + ^ E - ™ + V^) + 7E, (3.24) 

k=\ k=l ^ ^ k=l ^ ^ 

where the third, approximate, equality has been obtained by extending the upper limit of the 
sum from n to oo, which is correct up to terms of order C(e^"). The ensuing sum is evaluated 
in Appendix B under the assumption that r <S 1. Also, 7e = 0.577... is the Euler constant. 
Similarly for r > 2 we have (cf. Appendix B) 

n n -1 oo 1 /I iT^\^ 

k=l k=l \^ " ' fe=l v-^ " ' 
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where ({r) is the Riemann Zeta function. Note that the above summations are dominated by 
large values k = 0{n), as anticipated in the discussion preceding Eq. (3.16). Using Eqs. (3.24) 
and (3.25), one can do the summation over k in Eq. (3.23) to obtain 

l-a[fnPn{uj)] ^ -u;ln(l/T) -nln(l+Tu;) +lnr(l + a;). (3.26) 
Thus we finally arrive at 

^ ' J 27ri (1 + ™" ' ^ ^ 

c ^ ' 

which so far is valid for large n ^ 1/r but arbitrary Y . It is straightforward to check that 
this distribution satisfies indeed the large-n version of the master equation, i.e., the equation 
obtained from Eq. (2.3) after replacing /„ ~ l/r, as appropriate when l/r. 

Even though Eq. (3.27) is considerably simpler than the general form, it is still difficult to 
proceed without any further approximations, because of the presence of the Gamma function 
which has single poles on the negative real axis. Since we are interested in large values of n 
and Y , we can evaluate the integral in Eq. (3.27) using a saddle point approximation. As we 
shall see, the saddle point will occur at a small value of oj, so that we can replace the Gamma 
function by the lowest order terms in its expansion around cc; = 0:r(l+cj)~l — 7eCl^. (This 
means that only the sum Si in Eq. (3.24) and the t'^u terms of the sums in Eq. (3.25) are kept 
in the subsequent analysis.) We then obtain 

r da; e'^'^'' r 
Pn{Y) — ~ exp[F{uo)], (3.28) 

J 27ri (1 + ™)" ^27r|F"(u;o)| 

where we have defined the variable 

^^y-ln(l/r)-,._ ^3^^^ 

T 

and the function 

F{u) =ru;u -nln{l + rij), (3.30) 

while the contour integral over u has been evaluated in the saddle point approximation. For 
the latter to be justified, at least one of the conditions ur ^ 1 (i.e., Y Y^ = ln(l/r)) and 
riT 1 needs to be satisfied. Besides, the saddle point must obey |c<Jo| <^ 1, for consistency 
with the previous manipulations. Clearly, the saddle point occurs at 

n — v 

000^ , (3.31) 

SO the condition |c<Jo| ^ 1 is tantamount to \n — u\ <^ tu. Evaluating F and F" at the saddle 
point and substituting in Eq. (3.28) we finally obtain a Poisson distribution, as anticipated: 



Notice that the saddle point conditions require that the integration contour crosses perpendicularly the 
real axis. 
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Fig. 8. For high Y the probabiUties Pn{y) follow a Poisson distribution for large values of n. 



(3.32) 



As aforementioned, this approximation is valid so long as |n — z/| ^ tu. This range covers in- 
deed the "bulk" of the distribution at large Y (see also Fig. 8) : when summed over n within this 
range, Eq. (3.32) yields a total probability equal to 1 up to terms of order 0{e^'^^). Moreover, the 
lower limit nmin = f — tv = au h essentially the same as the upper limit (previously denoted 
as ricr) of the validity range of the approximation for P„ constructed in Sect. 3.1 (which, we 
recall, consists in preserving only the first term in the sum in Eq. (3.9)). Indeed, if one redefines 
ricr via the condition that these two approximations match with each other (up to prefactors) 
when n = n^, then one finds = (l — r + 0(r^))z/ ^ T^min- Thus, the two approximations 
for Pn at large Y that we have constructed in this paper are complementary to each other, in 
the sense that, together, they cover all the interesting values of n and they approximately match 
with each other at the borderline n^in = '^cr between their respective domains of validity in n. 

Now let us discuss some aspects of this distribution. The maximum of the distribution occurs 
at n* = z/ + 1/2 and the value at the maximum is P„* = Pmax = 1/ V^nu, while the width of 
the distribution is proportional to We need to say here that terms of order C(ln(l/r)) in 
the location of the maximum have not been kept, since this would requires a calculation of Si 
in Eq. (3.24) to 0{t) accuracy (see Appendix B). 

From Eq. (3.32) is is straightforward to deduce the average dipole number at high energy — 
one thus finds (n) = u, which is indeed the same as in Eq. (2.23) — and, more generally, all 
the A;-body dipole densities, which are readily obtained as 



when 



(3.33) 



This equation exhibits mean-field-like factorization, as expected for the bulk distribution at high 



^■^We should mention here that the calculation of {n^''^) is consistent with the validity range of the 
distribution (3.32) only so long as /c < 1/r. 
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Fig. 9. Factorization of expectation values holds only for quantities dominated by the bulk of the prob- 
ability distribution; here the ratio of the (normal ordered) dipole-pair number with respect to the dipole 
number squared which approaches unity for Y » Yc, in agreement with Eq. (3.33). 

density . Tliis is illustrated in Fig. 9. 

Let us now use the Poisson distribution (3.32) to compute the contribution of the bulk config- 
urations to the expectation value of the S'-matrix for the scattering of m external dipoles. One 
thus finds 



bulk 



fm 

bulk; 



(3.34) 



as expected from the mean field approximation, cf. Eq. (3.2). However, for m > 2, this bulk 
contribution is exponentially smaller (at large rz/ ~ F — Fc) than the respective contribution 
of the rare configurations with only few dipoles, as previously computed in Eq. (3.18). Also, 
even for m = 1, the result above is smaller by a large factor l/irv) <^ 1 than the previous 
result in Eq. (3.20), which, we remember, is due to the configurations with n ^ Ucr. Thus, the 
present calculation of the bulk distribution confirms our previous conclusion that the S'-matrix 
is dominated by the rare configurations which involve relatively few dipoles. 

The above results enable us to also compute the high-energy behaviour of Z(u, F) for u 
within the range a < u < 1, and thus complete our previous estimate, Eq. (3.19). Indeed, for 
large F and u > a, Z(u, F) is dominated by the typical configurations, distributed according 
to Eq. (3.32). By using the latter within Eq. (2.9), one easily finds 



Z{u, F) ^ exp 



1-u 



for F > Fc and a < m < 1. 



(3.35) 



The support of this function is strongly peaked near u = 1, within a distance \ —u ~ t/(Y — Y^) 
which becomes smaller and smaller with increasing energy and/or decreasing r. But if we let 



Notice that this is not the factorization that one finds in the dipole picture, i.e., in the absence of 
saturation; in that case one rather has (n^^)) ~ k\{n)^ at large Y [3]. 
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y ^ oo for fixed u < 1 (even arbitrarily close to 1), then Z{u,Y) vanishes exponentially, 
which is another manifestation of the 'black disk' limit in our model. 

3.3 More on boost invariance and the initial conditions 

Although the average S'-matrix is frame-independent, by construction, within the model un- 
der consideration, the physical picture of the collision — in the sense of relevant configurations 
— depends very much upon the choice of the frame, as we shall now demonstrate via some 
explicit calculations. (This dependence has been previously noticed in Refs. [1, 12].) To that 
aim, it is convenient to consider the (toy-model analog of) dipole-nucleus scattering; that is, 
we shall assume that, at F = 0, the target is made with A > 2 dipoles, whereas the projectile 
contains only one dipole: P^f (0) = 5nA and -Pm(O) = 5mi. Note that this is the same physical 
problem as considered in Sect. 3.1 — one onium starts as a single dipole in its rest frame, while 
the other one starts as a collection of exactly A dipoles (with A denoted as m in Sect. 3.1) — 
but the terminology that we use is now different (in the sense of interchanging what we call 
'target' and 'projectile'), since we feel that the present terminology is more natural in relation 
with a physical dipole-nucleus scattering. 

For the aforementioned initial conditions, we shall compute the average S'-matrix in two 
different frames: the rest frame of the target (nucleus) and, respectively, that of the projectile 
(dipole). For more clarity, we shall denote the results of these two calculations by (S*)^^^ and, 
respectively, (-S')^^^^. More general initial conditions will be briefly discussed towards the end. 

(a) Target (nucleus) rest frame. When the nucleus is at rest, the whole evolution is given to 
the dipole wavefunction. This is precisely the situation analyzed in Sect. 3.1, from which we 
simply quote here the final result (cf. Eq. (3.18) with m ^ A) : 

{S)i^A = ^ ^ A>2 (nucleus rest frame) . (3.36) 

As also discussed in Sect. 3.1, in this frame the average S'-matrix is dominated by rare config- 
urations in the evolved system (the 'dipole'). 

(b) Projectile (dipole) rest frame. This situation is new with respect to Sect. 3.1, in the sense 
that the whole evolution is now given to a system which starts with more than one dipole at 
y = (the 'nucleus'). For simplicity we shall consider first the case A — 2 and then return to 
generic values for A. From the arguments in Sect. 3.1, one may expect the average S'-matrix 
to be controlled by the rare configurations in the evolved nucleus which are dilute (but this 
expectation turns out to be wrong !), so let us first compute the respective contribution to {S). 
Following the steps given in Sect. 3.1, but now with a different initial condition, it is straight- 
forward to find that the distribution of the rare configurations reads 

Pt\y) ^ for 2<n<nc., (3.37) 

where is basically the same as before, i.e., n^r ~ [a /t)(Y — Y^. (Note, however, that the 
value of Yc depends upon A\ see the discussion around Eq. (3.40) below.) At a first glance, this 
result seems to imply that the average S-matrix will be proportional to c^^^^, which would 
be at variance with the previous result, Eq. (3.36), obtained in the nucleus rest frame and thus 
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would signal a violation of the boost invariance. However, this is not really the case, since the 
summation over the dilute configurations is not convergent anymore, rather it is dominated by its 
upper limit (which, we recall, is the borderline towards the bulk configurations). Specifically, 



Axl 



(3.38) 



n=2 



n=2 



where only the dominant exponential behaviour has been kept, and this comes from the terms 
with n ~ Tier, as anticipated. This should be contrasted with the calculation in the nucleus rest 
frame where the corresponding sum, cf. Eq. (3.17), is rapidly convergent. 

The above argument also shows that, in the present frame, (S) is rather dominated by the 
typical configurations in the evolved nucleus. Returning to the case of arbitrary A > 2, and 
following the procedure of Sect. 3.2 with the modified initial condition Pn(0) = SnA, one finds 
that the precise form of the Poisson distribution which describes the bulk is 



Pn{y) 



[n 



A)\ 



(3.39) 



with y = \Y — Yc{A)]/t, and where the critical rapidity Yc{A), determining the transition of 
the nucleus from the unsaturated "phase" to the saturated one, reads (for r <C 1) 



In ^ for At < 1 

At 



(3.40) 







for At > 1 . 



(Indeed, when At ^ 1, the nucleus is saturated already at F = 0. On the other hand, when 
At <C 1, the nucleus is dilute in the early stages of the evolution, when the average dipole 
number rises like (n) = Ae^ , cf. Eq. (2.23), until it reaches a saturation value of C(l/r); this 
happens when Y r^Yc with Yc as given above.) Given the probability distribution Eq. (3.39), it 
is straightforward to calculate the average ^'-matrix for the scattering of one dipole off the bulk 
distribution of the nucleus. One finds 



('S')axi 



At 



for At -C 1 
for At > 1 



(dipole rest frame) . (3.41) 



Thus is the final result in the dipole rest frame and is indeed in agreement with the corresponding 
one in the nucleus rest frame, as given in Eq. (3.36), so long as t <S 1 and A^l. 

Thus we see that, whether the (average) S'-matrix is dominated by rare configurations or the 
bulk distribution, really depends upon the frame that one decides to view the process in. This 
is not perhaps not so surprising, since it reflects the fact that the wavefunction of an evolved 
hadron is not a frame independent quantity. 

Notice also that there is some loss of accuracy when working with the bulk distribution rather 
than with the rare configurations: although the results (3.36) and (3.41) are consistent with each 
other within their validity ranges, the former is exact at high energy (i.e., it holds for any values 
of T and A), whereas the latter is only approximate (in our second calculation, based on the bulk 
distribution, we have not been able to obtain the expression of the prefactor when At ~ 1). 
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Now let us consider that the two systems are initially composed of A and B dipoles respec- 
tively, i.e. -P^(O) = and P^(0) = 5^3, and let us momentarily assume that A < B. From 
the previous discussion one understands that it is advantageous to work in the rest frame of 
the left mover. Then one can compute the asymptotic form of the average yS-matrix, since it is 
dominated by the rare configurations of the right mover wavefunction, and one finds 



AB 

e-^^^ for A<B 



(3.42) 

Ye-f^^ for A^B. 



The average S'-matrix for the special case A = B, has been obtained from the A-th equation of 
the hierarchy in Eq. (2.14), after {S)^^(^^_^i^ has been determined. Notice that Eq. (3.42) is an 
exact solution to this hierarchy. 

More generally, one could imagine that the initial wavefunctions are determined by some 
smooth distributions in the number of dipoles (rather than by a given number of dipoles). Then 
one can determine the S'-matrix, by averaging (S)^^^ over the initial conditions. So long as 
P^{0) and Pi(0) are non-zero, the asymptotic behaviour is simply given by the appropriate 
generalization of Eq. (3.20), namely 



(S) = E PA {0)Ph{0) (^)^xB ^ A^(0) A'(0) Ye-"- , (3.43) 

A,B=1 1 - <^ 

since the term A — B — 1 dominates the double sum when y — > oo. 



4 Correspondence between the toy model and high energy QCD 

Throughout this paper, we have emphasized similarities between the structure and the pre- 
dictions of the toy model and the known or expected properties of high-energy QCD. In what 
follows, we discuss this correspondence in a more systematic way. 

The factorization formula for the S'-matrix, Eq. (2.1), is reminiscent of the factorization 
schemes proposed within high-energy QCD in Refs. [3, 11,22,23,48], which have in com- 
mon to be symmetric between the projectile and the target. With the noticeable exception of 
the dipole factorization [3], which however fails to accommodate the saturation effects in the 
wavefunctions of the colliding hadrons, all the other schemes alluded to above are not written 
in terms of particle numbers, but rather in terms of gluon fields, sometimes represented (in the 
CGC formalism [37]) as classical color charges together with the color fields they radiate. It is 
therefore essential to establish the proper correspondence between gluons in QCD and particles 
in the toy model. 

This correspondence goes as follows: the 'dipoles' in the toy model correspond to s-channel 
gluons in QCD (so like the real color dipoles in Mueller's dipole picture), while r — 1 — a 
corresponds to gluons exchanged in the ^-channel. (The analog of r in QCD starts at order a^, 
corresponding to a two-gluon exchange between a pair of dipoles.) Hence, in an effective gluon 
language, the evolution described by Eqs. (2.3) and (2.5) contains vertices for gluon splitting 
in the s-channel — these are, of course, the vertices /„ — and also vertices for both splitting 
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and merging in the t-channel — each exchange r being sandwiched between two such vertices. 
Hence, this evolution constructs the analog of the 'Pomeron loops' in the onium wavefunction. 
Fully symmetric 'Pomeron loops' will develop, of course, in the ^-matrix for onium-onium 
scattering, as described by Eq. (2.1). 

With this identification, the toy-model hierarchy for the dipole scattering amplitudes, cf. 
Eqs. (2.14) or (2.19)-(2.21), corresponds to the QCD evolution equations with Pomeron loops, 
as constructed in Refs. [17-19]. To render this analogy more precise, let us focus first on the case 
where Pomeron loops are absent, and for which the toy model version is given in Eqs. (2. 17) and 
(2.18). We would like to argue that the continuum version of the master equation, Eq. (2.17), is 
the toy-model analog of the JIMWLK equation in QCD. The JIMWLK equation [7-9] is written 
for the color fields = a" radiated in the i-channel by the gluons produced in the s-channel 
by the high-energy evolution of the (right-moving) target. It is a second-order, functional, 
differential equation with respect to a and reads, schematically, 

where the transverse coordinates have not been shown. The quantity l^y[a] is the probability 
distribution for the fields (the analog of Pn{Y) of the toy model). The kernel x"*[q;] — 
the analog of the emission rate — is non-linear in a to all orders, via Wilson lines. These 
describe the multiple scattering of the gluon emitted in the s-channel in the last step of the 
evolution off the background field a (the t-channel gluons) produced by the previous steps of 
the evolution; this rescattering is similar to that included in /„ within the toy model. To complete 
the identification between Eqs. (4.1) and (2.17), one should recall that in the JIMWLK equation 
each s-channel gluon is allowed to radiate only two gluons in the i-channel. Hence, the second- 
order derivative w.r.t. a in Eq. (4.1) can be interpreted as a single derivative w.r.t. n (the number 
of gluons in the s-channel) and then Eq. (4.1) becomes indeed similar to Eq. (2.17). 

Strictly speaking, the JIMWLK equation does include some fluctuations, since it is Fokker- 
Planck equation for a random walk in the functional space of color fields [10]. However, these 
are merely color fluctuations which are suppressed at large Nc. The essential correlations as- 
sociated with gluon-number fluctuations are lost because of the impossibility to probe both 
(s-channel) gluons which are produced after one splitting. This is related to the fact that, as 
aforementioned, an s-channel gluon cannot radiate more than two t-channel fields, and thus 
cannot undergo multiple scattering off the projectile. By contrast, within the toy model, the 
dipoles are allowed to scatter multiply — both inside the target wavefunction and with the 
dipoles in the projectile — , hence the n-body correlations associated with splitting are fully 
taken into account. In the language of QCD, each s-channel gluon is allowed to absorb and ra- 
diate an arbitrary number of i-channel gluons, so like in the more general, 'self-dual', evolution 
described in Refs. [24,25]. 

In fact, as noticed in Ref. [33], the structure of the master equation for the toy model bears 
some formal resemblance to that of the 'diamond' Hamiltonian constructed in Refs. [24,25]. In 
particular, it shares with the latter the property of being self-dual, as it should, since the self- 
duality of the Hamiltonian is equivalent to the fact that the S'-matrix is boost-invariant, when 
the latter is given by a factorization formula like the one in Eq. (2.1) [22,23]. To see this, note 
that the master equation can be rewritten as (cf. Eqs. (2.3) and (2.5)) 
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dPn{Y) 




- yl — e 9" 




1-e- 



) P„(F) = H n\\na 



d_- 
dn. 



Pn{Y) , 



(4.2) 



dY 



where we have extended n to be a continuum variable and used g{n — 1) — exp{—d/dn)g{n) 
for a generic function g{n). As anticipated, the 'Hamiltonian' appearing in Eq. (4.2) is 'self- 
dual' [33], that is, it is invariant under the self-duality transformation which in the present 
context consists in exchanging 



and then reversing the order of the operators. Moreover, the presence of two types of exponen- 
tials — one involving the number of particles n and one involving the derivative d/dn — is 
reminiscent of the two types of Wilson lines which appear in the QCD Hamiltonian proposed in 
Refs. [24,25]. In QCD, n| In (j| is replaced by the color field produced by the s-channel gluons, 
and d/dnhy the (functional) derivative with respect to the color charge density of these gluons. 

Let us now return to the equations obeyed by the dipole scattering amplitudes and compare 
the structure of the fluctuation terms between the toy model and QCD. Consider first the equa- 
tion for (t^), as displayed in Eq. (2.20) for the toy model and, respectively, in Eqs. (2.7) and 
(2.8) of Ref. [18] for QCD. In both cases, the dominant fluctuation term is of the generic type 
r(t), i.e., it is linear in {t) and of order r ~ a^. Yet, the experience with the toy model tells us 
that the coefficient of the fluctuation term in QCD (cf. Eq. (2.7) in Ref. [18]) may actually be 
incomplete. Indeed, in the calculation of this term in Refs. [17-19] one has neglected the pos- 
sibility that the individual dipoles from the target undergo multiple scattering with the dipoles 
in the projectile; that is, in deriving the equation for (t^), the two external dipoles have been 
allowed to scatter only with two different dipoles from the target (but not also to scatter both 
off a same dipole). By contrast, in the toy model calculation, multiple scattering is included 
and is in fact responsible for an 0{1) contribution to the coefficient of fluctuation term r(t) in 
Eq. (2.20) : if that contribution was neglected, the ensuing coefficient would be twice as large 
than the correct one. More generally, if one considers the equation satisfied by (t^) within the 
toy model, then among the m — 1 (relevant) fluctuation terms which are included in this equation 
(cf. Sect. 2.2) — namely, T{t'^-'^), T^{t"'-'^), r'^-^it) — only the analog of the first term, 
r(t™^), has been so far included in the corresponding equation in QCD [18]. 

This discussion allows us to draw some lessons for QCD: most likely, the Pomeron loop 
equations in QCD at large N^. should be completed by including the effects of the multiple 
scattering of the individual target dipoles. This conclusion appears to be in conflict with some 
recent analyses within QCD [49, 50], which should be therefore carefully reexamined. Further- 
more, the more intricate structure for the fluctuation terms suggested by the toy model seems to 
prevent one from mapping the problem under study into a single Langevin equation. Note that, 
although the detailed structure of the fluctuation terms in QCD may be indeed more complicated 
than anticipated by the original analysis in Refs. [17-19,23], it is possible that the additional 
terms will not modify the qualitative asymptotic behaviour of the dipole scattering amplitudes 
at high energy and large Nc. As argued in Refs. [15, 17], this behaviour seems to be dictated 
simply by the stochastic nature of the evolution equations. 

An other interesting result of the toy model that we expect to extend to QCD as well is 
the fact that the average S'-matrix in the vicinity of the unitarity limit is dominated by rare 



—- < — > n\ In a 



(4.3) 
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fluctuations with only few gluons. This is in agreement with the arguments presented in that 
sense in Refs. [1,12], which also show that the exponential approach towards the black disk limit 
should be somehow faster in QCD, namely like (s) ~ exp(— F^), to be compared with the 
toy-model result (s) ~ exp(— F). (The additional factor of Y in QCD comes from the phase- 
space for diffusion in the dipole transverse sizes.) But what is most interesting about our present 
results is that the S'-matrix (s™) for the simultaneous scattering of several dipoles is even more 
strongly dominated by the rare fluctuations which are dilute, to the point that (s™) with m > 2 
is not at all sensitive to saturation effects in the target wavefunction, and (unlike (s)) it could 
have been simply computed within the dipole picture. It would be interesting to understand 
whether a similar feature holds in QCD as well. If so, this would mean that, for m > 2, (s"*) 
in high-energy QCD at large Nc can be reliably computed via numerical simulations within 
Mueller's dipole picture. 
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A Lack of boost invariance for the recombination process 

In this appendix, we shall consider more general (zero-dimensional) models in which in 
addition to particle splitting, one allows for recombination. The simplest model of this type 
is the reaction model A ^ AA, introduced in the context of QCD in Ref. [17], and further 
discussed in Refs. [32, 34, 35]. In that model, a particle can split into two (A — > AA) at a rate 
a per particle and, conversely, two particles can recombine into one (AA — > A) with a rate P 
per pair of particles. In what follows, we shall consider a more general version of this model, 
in which the rates for splitting (/„) and merging ((?„) are taken to be general functions of the 
number n of particles in the system. Hence, in such a process, particle number saturation can in 
principle occur via two different mechanisms: the saturation of the rate for particle splitting and 
the particle recombination. 

However, as we shall show in what follows, the recombination process cannot be made con- 
sistent with the boost invariance of the S'-matrix (within the factorization scheme of Eq. (2.1)) 
for any choice of the gn's. This finding, together with the fact that, within the context of the 
JIMWLK equation, gluon saturation occurs via the saturation of the emission rate, suggests that 
the models based on recombination may not closely resemble the QCD dynamics. This con- 
clusion is further supported by the observation in Ref. [26] that, if one attempts to interpret the 
Pomeron loop equations of QCD [17] as a reaction-diffusion process (in either the target or 
the projectile), then one is lead to introduce an effective 'recombination vertex' which has no 
definite sign, and hence no probabilistic interpretation. 

The master equation is obtained by adding to the r. h. s. of Eq. (2.3) the terms responsible 
for recombination. We thus obtain 



This is also the behaviour found by Salam, via numerical simulations of the onium-onium scattering 
within the context of the dipole picture [1,4]. 
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dP (Y) 

-^r^ = fn-l Pn^liY) - fnPn{Y) + Qn+l Pn+l{Y) - gnPn{Y), (A.l) 

where the functions and gn should be chosen so as to ensure the boost invariance of the 
^-matrix and satisfy the boundary conditions fi = a and gi = 0. By requiring the average 
S'-matrix in Eq. (2.1) to be independent of Yq, we easily obtain the following constraint 

U{1 - a") - fn{l - a^) - gmi'J-'' - 1) + gni'J-"' - 1) = 0. (A.2) 

A priori, there are two different ways in which this constraint could be satisfied: (i) the inde- 
pendent cancellation of the /-terms and -terms, and (ii) the mutual cancellation between the 
two types of terms. Both are excluded, as we now show. 

(i) By requiring the independent cancellation of the /-terms and the g'-terms in Eq. (A.2), 
one finds 

/„ = c(l-(7") and g„ = d(a-"-l), (A.3) 

where the constants c and d are determined by the boundary conditions at n = 1 . For /„ we 
recover Eq. (2.5), but the condition gi = can be satisfied only by choosing d = 0, and hence 
Qn — O for any n. 

(ii) In order for the /-terms to cancel against the g'-terms in Eq. (A.2) one needs ^® 

/„ = c(a-"-l) and 5„ = -c(l-a"). (A.4) 

Depending on the sign of the constant c, either the /„'s or the (?„'s will be negative, and thus there 
is no probabilistic interpretation — a situation reminiscent of the difficulty met in Ref. [26]. 
The standard reaction-diffusion process A ^ AA is defined by (see, e.g., [17]) 

fn = an and gn = Pn{n - l)/2, (A.5) 

and thus it does not belong to the case (ii); this process has a well-defined probabilistic 
interpretation, but is inconsistent with the boost invariance of the ^'-matrix in the factorization 
scheme of Eq. (2.1), which includes multiple scattering for the individual particles . 

To better appreciate this difficulty, let us explicitly compute, within the reaction model of 
Eq. (A.5), the high-energy limit of the average S-matrix for 'dipole-nucleus scattering' — 
i.e., for the scattering between two onia which in their respective rest frames reduce to a single 
dipole in the case of the projectile and, respectively, to a set of A dipoles, with A > 1, in the 
case of the target — by separately working in the target rest frame and in the projectile rest 
frame. We anticipate that the final results will be non-zero in both cases (in agreement with the 
'grey disk limit' reported in Ref. [35]), but different for the two calculations. Only one of these 



^^From the perspective of Eq. (2.1), this would correspond to a scenario in which the emission of a 
dipole in one onium can be reinterpreted, after a shift dY in Yq,2l?,sl recombination in the other onium. 

In Ref. [23], the Pomeron loop equations in QCD at large- A^^ have been given an effective interpreta- 
tion in terms of a reaction-diffusion process with 'BFKL Pomerons', which appears to be self-dual and 
thus formally consistent with boost invariance. Note, however, that in that construction, one has given up 
multiple scattering, that is, one has used a different factorization scheme in which the individual dipoles 
(from the target and the projectile) are allowed to scatter only once. 
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results will correspond to that found in Ref. [35], which has adopted the point of view of the 
projectile evolution. 

Before we proceed, let us notice that, within the reaction model of Eq. (A.5), the quantities 
a, P and cr are a priori independent parameters (the first two of them referring to the evolution 
of the onium wavefunction, and the third one characterizing the elementary dipole-dipole scat- 
tering), but in view of the correspondence with QCD one needs to assume that the ratio (3 /a 
is of the same order as r = 1 — a, namely they are both of 0{a'^) (in particular, this implies 
/? <^ a). To facilitate the comparison between this model and the one studied in this paper, we 
denote /3/2a = tq, with tq ~ r <i; 1. 

We start by summarizing the formulae giving {S)y in the considered frames: 

(a) Target (nucleus) rest frame (i.e., projectile evolution). This is the case Yq — Y in 
Eq. (2.1), which then implies: 

(5)^ = j:P^iY) a^^ = Z'^iu = a^,Y), (A.6) 

m 

where P^(0) = dmi, so that Z^{u, Y = 0) = u. 

(b) Projectile (dipole) rest frame (i.e., target evolution). In this case, Yq — and Eq. (2.1) 
implies: 

{S)y = E Pniy) = Z^iu ^C7,Y), (A.7) 

n 

where this time -P^f (0) = 6nA, implying Z^ {u, Y = 0) = u^. 

We have related here {S)y to the generating functional Z{u,Y), cf. Eq. (2.13), since as 
shown in Ref. [35] the latter is the quantity which is most conveniently evaluated in the context 
of the reaction model. Specifically, it is first straightforward to use Eq. (A.l) with the vertices 
/„ and Qn from Eq. (A.5) to deduce the following evolution equation for Z{u, Y) in the reaction 
model (to be compared to Eq. (2. 11)) : 

, ,{ dZ (3 d^Z] 
_=„(1_„)|_„_ + |_|. (A.8) 

This equation is the starting point of the analysis in Ref. [35]. As shown there, and it can be also 
easily verified by inspection, this equation admits the following, non-trivial, fixed point 

= e 1 1 ' (^0 = /^/2«)' (A.9) 

which, as anticipated by our notations, is the same as the high-energy limit of the physical 
generating functional (i.e., the solution to Eq. (A.8) with physical initial conditions) : Z{u, Y — > 
oo) Z^{u) for any initial conditions. (This is demonstrated in Ref. [35].) This should be 

contrasted to the corresponding behaviour of the model studied by us here, where we have seen 
that Z(u, Y) vanishes exponentially when F ^ oo for any u < 1 (cf. Eqs. (3.19) and (3.35)). 

By substituting Eq. (A. 9) in Eqs. (A.6) and (A.7), it is easy to compute the asymptotic values 
of the (average) ^'-matrix in the two Lorentz frames under consideration. Explicitly, one obtains 
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{S)l = e-(-^^)/- - e-V- , (A.10) 

where the second, approximate, equality holds since tq <^ 1. The quantity 1 — a"^ can be 
recognized as the original value of the scattering amplitude at y = 0: (T')q — 1 — {S)^ — l—a"^. 
It is not hard to show that {S)^ < {S)^ — for any positive value of tq and any value of 
such that < < 1. Thus, in this case, one reaches the natural conclusion that the value of 
the ^-matrix becomes smaller at high energy. 
On the other hand, 

(t/to _ 1 

= — ~ e-^/^° - e-^/^°, (A.ll) 

where the exponent t/tq is of order one. Clearly, for any A > 1, the two above results are 
different from each other, namely, one has 1 - CT^ > 1 - a = r and hence (5)^ < (S)^. 

Moreover, when A is large enough, the results above predict that the asymptotic S'-matrix 
at high energy may become larger than the original .S-matrix {S)^ = at F = 0. Indeed, by 
using I In cr I ~ r, one can easily check that (S)^ > as soon as ^4 > I/tq. 



B A useful sum 

In this Appendix we calculate the following sum (with a — 1 — r) 



oo 



«5i(r)^rEY3^. (B.l) 

n=l ^ " 



Up to the T prefactor, this is equal to the double sum J2m,n=i which cannot be written in 
terms of known functions. Nevertheless, we can obtain an analytic expression for r ^ 1, by 
keeping the first few terms in a series expansion. Here we shall work up to order 0{t). (Of 
course, the summation can be easily evaluated numerically, since it is rapidly convergent for 
any a < 1.) 

One can convert the sum into an integral by making use of the Euler-McLaurin summation 
formula, which may be written as 

Jdnh{n) + 1 Ml) + 2 E ^"^ ^ h^''-'\l). (B.2) 

n=i ■{ ^ k=i y^^) 

With h{n) — TO-"/ (1 — cr"), it is straightforward to calculate the integral and obtain 

oo 

/ dn = — - = ln rln- + Or2. (B.3) 

y l-cr" In(l-T) T 2 T ^ ' 

The integral contains the dominant contribution for small r, which is equal to ln(l/T). It is 
trivial to obtain the contribution of the second term in Eq. (B.2) which is 
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Ml) 



1 r 

2 ~ 2' 



(B.4) 



For the contribution of the third term in Eq. (B.2) we need the (2fc — l)-th derivative of h{n) at 
n = 1 which, to the order of interest, reads 



/i(2fc-i)(l) = -T{2k) (1 - - ) + 0(r 



(B.5) 



Notice that we have to sum a series whose terms are growing factorially. However, the terms 
have alternating signs and the series is Borel summable. Thus we can write the third term as 



-2 1 



E 



C(2A;)r(2A:) 



k=l 

oo 



(27r)2fc 



2k 



- coth - 



7E-2 



r 

1 - - 
2 



(B.6) 



Adding the three contributions obtained in Eqs. (B.3), (B.4) and (B.6) we finally arrive at 

5i(r) = In i + 7E - ^ r In - - r + ^(r^). (B.7) 

Similarly one can obtain the other sums encountered in Eq. (3.25). One finds 

SJt) = T — ^ = C(^) + 0{t) for r > 2. (B.8) 



n- 
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